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SUMMARY 


OBJECTIVE 

Fit  irregularly  spaced  data  with  continuous  mathematically  differentiable  hmcdons. 

RESULTS 

A  numerical  process  was  developed  accomplishing  the  objective.  Applications  to  mapping  of 
ionospheric  parameters  and  HF  ray  tracing  in  an  empirical  ionosphere  were  successfuUy  tested. 


RECOMMENDATIONS 

Since  the  functions  used  in  the  data-fitting  procedure  are  nonlinear,  more  testing  with  a  variety  of 
data  sets  may  be  required  for  optimalization  and  error  estimates. 
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INTRODUCTION 


Frequently,  data  are  collected  at  irregularly  spaced  grid  points,  and  for  a  variety  of  reasons,  there 
is  usually  a  need  for  interpolation;  for  example,  a  two-dimensional  graphical  representation  in  the 
form  of  maps  showing  contour  lines  of  constant  values  (isobars  in  weather  maps).  In  many  cases,  this 
is  not  a  trivial  task,  especially  when  the  grid  points  are  irregular  and  large  data  gaps  exist.  This  prob¬ 
lem,  as  it  exists  for  ionospheric  data,  was  discussed  by  Rush  et  al.  (1983).  If  the  data  can  be  fitted 
with  a  continuous  function  that  has  continuous  partial  derivatives,  the  computation  of  contour  lines  is 
relatively  easy. 

The  electron  density  in  the  ionosphere  at  a  given  time  is  a  function  of  latitude,  longitude,  and 
height.  Horizontal  gradients  of  the  electron  density  can  cause  significant  deviation  of  the  ray  path  dis¬ 
tance  and  direction  from  that  expected  under  the  assumption  of  horizontal  stratification  of  the  iono¬ 
sphere.  To  improve  propagation  predictions  under  such  circumstances,  ray  tracing  studies  that  require 
the  knowledge  of  the  density  and  its  gradient  for  every  point  along  each  ray  path  are  needed.  Similar 
remarks  apply  to  the  vertical  sovmding  of  the  ionosphere,  where  horizontal  gradients  cause  deviation 
from  vertical  propagation  and  limit  the  accuracy  of  electron  density  profile  calculations  unless  angle-of- 
arrival  data  are  available  and  corresponding  corrections  can  be  obtained  by  means  of  ray  tracing 
(Paul,  1985;  Paul,  1989;  Wright,  1990). 

In  the  following,  we  introduce  a  process  for  fitting  discrete  data  given  at  irregular  grid  points  with 
one  function  that  is  continuous  and  has  finite  continuous  first-  and  higher-order  derivatives  every¬ 
where. 

THE  FTITING  FUNCTION 

The  fitting  function  process  is  based  on  a  simple  concept.  Any  function  with  the  following  proper¬ 
ties  can  be  used: 

a.  The  one  or  multidimensional  function  gi  is  of  unit  value  at  a  given  location  indicated  by  the 
index  i;  it  decreases  and  approaches  zero  value  asymptotically  with  increasing  distance  from 
this  location  (e.g.,  Gaussian). 

b.  The  rate  of  the  isotropic  or  anisotropic  decrease  of  gi  is  controlled  by  one  or  several  control 
distances. 

c.  At  least  the  first-order  partial  derivatives  of  gi  exist,  are  continuous,  and  have  finite  value. 

For  each  location  where  data  are  given,  a  “bumplet"  is  defined  by 


//=  1+0/  •  gi. 


(1) 


The  functional  description  of  the  data  is  then  given  by  the  product  of  one  normalizing  factor  and 
all  the  bumplets  fi  (one  for  each  data  point).  For  the  two-dimensional  case,  we  then  have 


P(x,y)  =  Pon/i(x,y). 


(2) 


The  parameter  oi  in  each  term  is  obtained  by  an  iteration  process  to  fit  the  given  data. 
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One  type  of  function  that  meets  the  requirements  specified  above  (Paul.  Smith,  and  Wright,  1968) 
is  given  by 


(3) 


where 


Xri 


x-xi 

rxl 


yn 


y-yt 

ryt 


(4) 


and  m  is  an  even  integer  number.  The  control  distances  in  the  x  and  y  direction,  re^ctively,  are 
fxi  and  Tyi.  For  a  curve  given  by 


x^r  +  yS  =  1 


(5) 


(an  ellipse  for  m=2),  the  value  of  the  function  gi  is  0.5.  The  shape  of  this  curve  depends  on  the 
exponent  m,  but  it  intersects  the  x-  and  y-axis  at  the  same  points  given  by  the  corresponding  control 
distances  for  all  even  values  of  m.  Figure  1  shows  an  example  of  the  function  g  (one-dimensional 
case)  for  two  values  of  the  exponents  m-2  and  m=4.  Equation  (3)  can  be  modified  to  permit  the  use 
of  rational  or  irrational  numbers  for  the  exponent  m,  if  the  absolute  values  of  the  relative  distances 
Xri  and  yri  are  used.  In  the  remainder  of  this  report,  it  will  be  assumed  that 


m  =  2 


(6) 


unless  otherwise  specified.  Generalization  to  the  n-dimensional  case  and  higher  values  of  m  is  straight¬ 
forward. 

After  the  a/  parameters  have  been  determined  in  an  iteration  process  (to  be  discussed  later),  the 
partial  derivatives  can  be  obtained  to 


^  -  _2-  /<  -  ^ 

6x  ~  fxi  I  +  x%  +  y%’ 

^  -  _2-  ^  ^ 

6y  ~  ~  ryi  1  +  +  yj,  ■ 

The  partial  derivatives  of  the  function  p  can  then  be  e;q>ressed  in  the  forms 

^  -D  •  Y  1  ^ 

6x  ^  fi  6x’ 


^  =D  •  \  ^ 

dy  ^  fi  6y' 


(7) 

(8) 

(9) 

(10) 
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THE  ITERATION  PROCESS  FOR  DATA  FITTING 


In  the  two-dimensional  case,  the  function  p  defined  in  equation  (2)  may  contain  as  many  as  3n-t-l 
unknown  parameters  if  n  data  points  are  avaUable.  This  indicates  that  the  process  is  flexible  and  at 
the  same  time  not  unique.  The  user  wUl  have  to  select  some  of  the  parameters  on  the  basis  of  what 
appears  reasonable.  In  most  cases,  it  is  reasonable  to  set  the  control  distances  and  ryi  equal  to  or 
proportional  to  the  distance  between  the  reference  point  and  its  nearest  neighbor. 

The  choice  of  po  is  not  critical,  and  in  many  cases,  the  average  value  of  the  data  can  be  used.  If 
the  data  set  should  contain  positive  and  negative  values,  it  may  be  better  to  add  a  constant  value  to 
the  data  and  subtract  this  constant  from  the  expression  on  the  right-hand  side  of  equation  (2).  By 
using  this  procedure,  a  zero  value  of  Po,  if  taken  as  the  mean  value  of  the  data,  can  be  avoided. 
Equation  (2)  also  shows  that  the  fitting  function  approaches  the  value  po  in  areas  far  away  from  loca¬ 
tions  where  data  had  been  recorded. 

This  now  leaves  the  n  parameters  a/  to  be  determined  from  the  data  A'.  In  a  first  approximation 
for  the  iteration  process,  we  set 

1 

(„) 


We  then  compute  values  pi  for  the  grid  points  by  using  equation  (2).  From  the  difference 
between  those  values  and  the  data,  a  correction  of  doi  can  be  derived  by 


doi 


Dj  -  Pi 
Po 


(12) 


Usually,  only  a  fraction  of  this  value  is  applied  as  a  correction  to  the  a/,  and  the  iteration  process 
is  repeated.  In  this  way,  oscillatory  behavior  and  poor  convergence  can  be  avoided.  Depending  on  the 
accuracy  requirements  and  the  variation  and  distribution  of  the  data,  a  small  number  of  iterations  usu¬ 
ally  yields  a  good  agreement  between  the  fitting  function  and  the  input  data.  Since  the  fitting  function 
is  not  linear,  it  is  not  possible  to  make  more  precise  statements  about  the  convergence  of  the  iteration 
process. 

Figure  2  shows  an  example  of  connecting  six  data  points,  given  at  irregular  intervals,  with  a 
smooth  curve  p(x)  by  using  this  process.  Also  shown  in  figure  2  is  the  first  derivative  dpidx  of  the  fit¬ 
ting  function. 

One  advantage  of  this  fitting  procedure  is  its  flexibility.  That  is,  local  adjustments  do  not  cause 
large  global  changes,  as  demonstrated  in  figure  3.  The  dashed  curve  is  the  same  as  that  in  figure  2. 
For  the  solid  curve,  the  same  data  were  used  with  the  exception  that  the  third  data  point  was  changed 
from  7.5  to  2.S.  Significant  changes  are  visible  in  the  range  between  2  and  S,  while  there  is  very  little 
change  in  the  remainder  of  the  fitting  curve. 
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APPLICATION  TO  IONOSPHERIC  DATA 


As  mentioned  earlier,  a  fitting  process  as  described  here  can  be  used  in  several  different  applica¬ 
tions  with  ionospheric  data.  For  example,  for  several  observed  parameters,  a  worldwide  description  in 
the  form  of  maps  is  desirable.  Furthermore,  ray  tracing  studies  require  continuity  of  the  electron  den¬ 
sity  and  its  gradients  in  three  dimensions.  Table  1  lists  the  monthly  median  of  the  ordinary  critical 
frequency  of  the  F-region  foF2  for  December  1960  at  0700  hours  UT  for  a  number  of  stations  in 
Europe.  This  area  and  the  particular  period  were  selected  for  two  reasons: 

1.  Europe  had  the  highest  density  of  observatories. 

2.  On  the  average,  the  strongest  electron  density  gradients  are  observed  in  the  sunrise  area  in 
winter  time  when  solar  activity  is  high. 

A  flat-earth  approximation  is  used  in  this  example  by  projecting  the  location  of  the  observatories 
into  a  plane  tangent  to  the  surface  of  the  earth  at  a  latitude  of  SO  degrees  and  a  longitude  of  20 
degrees.  The  corresponding  Cartesian  coordinates  in  this  plane  are  listed  in  the  last  two  columns  of 
table  1.  First,  a  map  of  foF2  was  generated  as  shown  in  the  form  of  contour  lines  of  constant  foF2 
in  figure  4  (solid  curves).  The  locations  of  the  stations  are  marked  by  asterisks.  The  contour  lines  in 
the  upper  left  part  of  the  hgure  show  the  expected  orientation  parallel  to  the  daylight  boundary.  This 
is  not  true  for  the  lower  right  part  of  this  figure,  where  virtually  no  data  were  recorded  and  the  fitting 
function  approaches  the  average  value  po-  A  second  map  was  generated  by  adding  one  artificial  data 
point  as  listed  last  in  table  1.  The  resulting  contour  lines  (dashed  lines)  are  quite  different  in  the 
vicinity  of  this  additional  point,  but  remain  virtually  unchanged  in  the  area  of  high  data  density.  This 
example  again  shows  the  flexibility  of  this  method  for  local  corrections  or  updating.  Figure  S  shows 
foF2  contour  lines  for  0800  hours,  1  hour  later  than  figure  4.  The  orientation  of  the  contour  lines 
changed  little,  but  their  values  increased  in  time  as  expected.  The  shallow  maximum  in  the  lower  cen¬ 
ter  (the  dashed  lines  correspond  to  a  maximum  plasma  frequency  of  8.5  MHz)  is  caused  by  the  lack 
of  data. 

Table  2  gives  some  details  of  the  iteration  process  for  producing  this  map.  The  first  colunm  lists 
the  foF2  values  and  the  next  columns  show  the  amplitude  and  the  differences  between  computed  and 
input  data  after  a  number  of  iterations.  The  last  two  columns  give  the  station  coordinates.  The  largest 
deviation  after  40  iterations  is  7  kHz,  which  for  all  practical  purposes,  should  meet  every  accuracy 
requirement.  After  40  iterations,  a  few  of  the  amplitudes  become  small,  as  can  be  seen  in  lines  S  and 
6  and  especially  in  line  17  of  table  2.  This  means  that  almost  the  same  fitting  function  can  be  pro¬ 
duced  without  the  three  input  data  listed  in  those  lines.  This  could  result  in  significant  savings  where 
computer  time  is  a  critical  issue. 

If  corresponding  maps  for  hm,  the  height  of  the  maximum  of  the  F-layer,  and  ym,  its  half-thick¬ 
ness,  could  be  produced,  a  parabolic  model  of  the  F-layer  suitable  for  ray  tracing  could  be  con¬ 
structed  since  the  plasma  frequency  and  its  gradient  would  be  defined  everywhere  in  the  spatial 
domain  of  interest.  Unfortunately,  many  fewer  hm  and  ym  data  than  foF2s  are  available  since  their 
production  requires  some  form  of  electron  density  profile  calculation,  a  much  more  complicated 
process  than  the  scaling  of  foF2  from  ionograms. 

An  alternative  approach  is  the  use  of  observed  or  theoretical  data  in  combination  with  reasonable 
assumptions  for  the  construction  of  an  ionospheric  model  suitable  for  ray  tracing.  Such  a  process  is 
outlined  in  the  following  paragraphs. 

Figure  6  shows  sue  profiles  of  plasma  frequency  as  a  function  of  height.  These  can  be  interpreted 
as  a  temporal  sequence  during  the  sunrise  period  at  a  fixed  location  or  sample  profiles  recorded  at  the 
same  time  at  different  locations  in  the  sunrise  zone.  Three  of  the  profiles  show  only  an  F-layer,  with 
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foF2  increasing  and  hm  decreasing  from  left  to  right.  The  other  three  profiles  include  an  E-layer  and 
a  minimum  (valley)  between  the  E-  and  the  F-layer.  Each  profile  is  defined  by  a  few  data  points 
(marked  by  Xs).  The  data  are  again  connected  by  the  type  of  fitting  function  given  in  equation  (2), 
which  also  guarantees  the  continuity  of  the  slope  at  each  point  of  the  profile.  This  facilitates  the  com¬ 
putation  of  the  virtual  heights  for  each  profile. 

A  truly  three-dimensional  ionospheric  model  was  generated  from  seven  profiles  similar  to  those 
shown  in  figure  6.  This  model  is  defined  by  34  data  points.  The  plasma  frequencies  and  their  coordi¬ 
nates  are  listed  in  table  3.  This  data  set  shows  that  typical  horizontal  scale  sizes  are  much  larger  than 
vertical  scale  sizes.  Therefore,  a  ratio  of  4:1  was  used  for  the  ratio  of  the  horizontal  to  the  vertical 
control  distances.  Each  of  the  seven  profiles  shows  a  negative  plasma  frequency  at  the  lower  end  of 
the  height  range.  This  is  done  to  ensure  that  the  fitting  function  reaches  a  zero  value  at  the  lower 
boundary  of  the  model  since  this  surface  defines  the  lower  edge  of  the  ionosphere.  The  fitting  func¬ 
tion  is  defined  everywhere,  but  only  that  portion  between  the  lower  boundary  where  the  plasma  fre¬ 
quency  is  zero  and  some  upper  boundary,  usually  the  height  of  maximum  plasma  frequency,  is  used. 

It  is  assumed,  at  least  for  ray  uacing  purposes,  that  the  area  between  the  ground  and  the  lower 
boundary  is  free  space  with  zero  plasma  frequency. 

The  complexity  of  such  a  model  ionosphere  is  demonstrated  in  figures  7  and  8  by  partial  vertical 
cross  sections  in  the  xz  plane  at  y  =  100  km  and  in  the  yz  plane  at  x  =  ISO  km,  respectively.  The 
solid  lines  are  contour  lines  of  constant  plasma  frequency  (in  MHz).  The  plasma  frequencies  indicated 
by  the  dashed  lines  differ  from  the  values  of  the  neighboring  solid  lines  by  0.5  MHz. 

Figures  9,  10,  and  11  show  some  ray  tracing  results  (earth’s  magnetic  field  neglected)  through 
such  an  ionosphere.  It  is  not  surprising  that,  under  those  conditions,  the  ray  path  is  strongly  asymmet¬ 
ric.  The  projection  into  the  ground  plane  (figure  11)  indicates  a  large  azimuthal  deviation  between  the 
angle  of  arrival  at  the  receiver  site  and  the  true  direction  to  the  transmitter.  It  also  shows  that  this 
azimuth  error  can  become  quite  different  when  receiver  and  transmitter  locations  are  interchanged.  A 
homing  procedure,  where  the  up-leg  and  down-leg  of  the  path  become  identical,  was  used  to  obtain 
the  ray  path  shown  in  figure  12.  As  expected,  the  deviation  from  vertical  propagation  is  rather  large. 
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DISCUSSION  AND  CONCLUSIONS 


As  mentioned  earlier,  the  fitting  process  is  not  unique  but  offers  high  flexibility.  This  will  be  Ulus- 
trated  by  one  more  example.  The  same  data  were  used  as  in  figure  2,  but  the  value  of  th  '  normaliz¬ 
ing  factor  po  was  set  equal  to  10.0.  The  results  for  two  different  sets  of  parameters  are  she  vn  in 
figure  13.  In  one  case  (solid  lines),  a  high  value  for  the  exponent  in  equation  (3)  was  used  (m  =  16) 
together  with  a  short  control  radius  (0.3  times  the  distance  to  the  closest  point).  In  the  other  case 
(dashed  lines),  the  exponent  was  low  (m  =  2),  and  the  entire  distance  to  the  nearest  point  (r  =  1.0) 
was  used  as  a  control  radius.  Either  case  could  be  representative  for  a  physical  process.  For  example, 
assume  that  during  daytime  the  output  of  a  light  meter  is  recorded  while  individual,  well-separated 
clouds  of  different  sizes  and  densities  move  across  the  sky.  Between  clouds,  the  brightness  is  constant 
(represented  by  a  value  of  10.0  in  figure  13),  then  drops  during  a  short  transition  time  to  a  lower 
value  (depending  on  the  size  and  the  density  of  the  cloud)  and  remains  approximately  constant  again 
until  the  cloud  moves  out.  The  second  (dashed)  curve  shown  in  figure  13  may  be  representative  of, 
for  example,  the  tidal  changes  of  the  sea  level.  Both  curves  fit  the  same  six  given  data  points.  Know¬ 
ing  the  physical  processes  involved,  it  is  obvious  that  the  dashed  curve  would  be  a  poor  approximation 
to  the  change  of  brightness  in  the  example  given  above,  and  it  would  be  absolutely  unrealistic  to  use 
the  solid  curve  for  a  description  of  the  tidal  variation.  It  can  be  argued  that,  implicitly,  there  is  addi¬ 
tional  information  available  for  the  solid  curve,  namely  the  background  value  for  full  brightness  po- 
The  main  difference  between  the  two  curves,  however,  results  from  the  fact  that,  in  the  case  repre¬ 
sented  by  the  solid  curve,  the  function  p  approaches  the  value  of  its  normalizing  factor  po  quickly  at 
distances  beyond  the  control  distance  because  the  control  distances  were  chosen  to  be  rather  short 
compared  with  the  distance  to  the  nearest  neighboring  point.  For  the  dashed  curve,  however,  when 
the  abscissa  moves  beyond  the  control  distance  of  one  datum,  it  is  already  within  (or  at  least  close  to) 
the  control  distance  of  the  next  data  point  and  carmot  reach  the  background  value,  unless  the  next 
datum  is  equal  or  close  to  po. 

The  fitting  process  described  here  features  weak  coupling  from  one  interval  to  the  next  in  contrast 
to  fits  using  spline  function,  where  the  coupling  is  strong.  Similar  remarks  apply  to  fitting  with 
sinusoids  or  polynomials  because  local  corrections  or  adjustments  of  the  data  cause  changes  over  the 
entire  fitting  range.  Looking  at  geophysical  data,  we  find  the  entire  spectrum  of  situations,  from  one 
driving  force  controlling  the  worldwide  response  such  as  the  changes  of  the  gravitational  field  causing 
ocean,  atmospheric,  and  solid  earth  tides,  to  very  local  events  that  leave  the  remainder  of  the  globe 
unchanged;  e.g.,  tornadoes.  In  between,  we  have  relatively  large  but  still  local  features;  e.g.,  low  or 
high  barometric  pressure  systems,  areas  of  rainfall,  and  magnetic  storm  effects  in  the  ionosphere.  The 
functional  variation  of  the  solid  curve  in  figure  13  may  also  be  representative  for  the  variation  of 
mechanical  or  electrical  properties  of  layered  rock  formations.  Depending  on  the  type  of  data,  a  local 
modeling  or  mapping  may  be  more  adequate  and  more  economical  than  an  attempt  at  global  repre¬ 
sentation. 

It  should  be  mentioned  that  the  weak  coupling  of  the  fitting  process  outlined  here  also  means  that 
this  procedure  may  not  be  suitable  for  extrapolation  or  prediction.  This  can  be  seen  in  figure  14 
where  the  data  set  used  in  figure  2  and  figure  13  for  m  =  2  are  compared.  The  difference  here  is 
that,  for  the  solid  curve,  the  average  value  (4.75)  of  the  data  was  used  for  po,  while  the  dashed  curve 
is  identical  to  the  corresponding  curve  of  figure  13  with  po  =  10.0.  There  is  a  significant  difference 
between  the  two  curves  to  the  left  side  of  the  lowest  data  point  and  a  small  difference  between  the 
fourth  and  the  fifth  data  points.  As  mentioned  above,  the  function  approaches  the  value  po  for  val¬ 
ues  far  away  from  data;  that  explains  the  different  behavior  of  the  two  curves  on  the  left  side  of  the 
figure  and  also  their  difference  between  the  fourth  and  the  fifth  data  points,  where  there  is  a  relatively 
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large  gap.  This  example  also  demonstrates  that  the  Nyquist  theorem,  which  requires  a  minimum  sam¬ 
pling  rate  if  aliasing  is  to  be  avoided,  does  not  guarantee  a  tmique  interpolation  between  consecutive 
sampled  values,  unless  it  is  known  that  the  sampled  function  is  a  superposition  of  sinusoids. 
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Figure  2.  Fit  of  six  irregularly  spaced  points  with  the  function  p(x)  and 
the  first  derivative  dp(x)ldx. 
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Figure  4.  Contour  lines  of  constant  maximum  plasma  frequency  foF2  (solid 
lines)  derived  from  dau  observed  at  locations  indicated  by  stars.  The  figure 
also  shows  the  effect  of  adding  one  artihcial  data  point  (dashed  lines). 


Median  foF2,  08.00  UT,  Dec.  1960,  Europe 


Figure  5.  Contour  lines  of  constant  maximum  plasma  frequency  1  hour  later 
than  those  shown  in  figure  4. 


SEQUENCE  OF  VERTICAL  PROFILES  FOR  SUNRISE  MODEL 


Figure  6.  Plasma  frequency  profiles  for  sunrise  conditions. 
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CONTOUR  LINES  OF  CONSTANT  PLASMA  FREQUENCY 


Figure  7.  Vertical  cross  section  (east-west  plane)  through  the 
lower  part  of  a  model  ionosphere  in  the  sunrise  zone. 


YZ  PLANE.  X  =  t50  km 


Figure  8.  Vertical  cross  section  (north-south  plane)  through  the  lower 
part  of  a  model  ionosphere  in  the  sunrise  zone. 


XY  PLANE. 


PR  =  4.00  MHz 


Figure  11.  Ground  projection  of  a  ray  path  through  the  sunrise 
model  ionosphere. 


X2  PLANE,  PR  *4.00  MHz 


Figure  12.  Deviation  of  a  ray  path  from  vertical  propagation. 
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M  =  16/2  R  =  0.3/1  (. 


Figure  13.  Fitting  a  set  of  data  with  the  same  type  of  functions  with 
different  parameters  (control  distance  and  exponent). 


M  =  2  R=1.0  PO  =  20.0/4.75 


Figure  14.  Effect  of  changing  the  normalization  factor. 
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Table  1.  Monthly  median  foF2,  December  1960,  0700  UT. 


Station 

Lat. 

Long. 

foF2 

X 

Y 

Inverness 

57.4 

355.8 

2.2 

-1463.5 

1093.8 

Juliusruh 

54.6 

13.4 

4.5 

-426.5 

532.6 

De  Bilt 

52.1 

5.2 

3.4 

-1013.5 

337.5 

Slough 

51.4 

359.4 

2.6 

-1435.4 

359.6 

Lindau 

51.6 

10.1 

4.5 

-684.6 

224.4 

Pruhonice 

50.0 

14.6 

5.6 

-386.0 

13.9 

Dourbes 

50.8 

4.3 

3.3 

-1106.3 

207.1 

Saclay 

48.1 

2.3 

3.4 

-1320.9 

-58.1 

Freiburg 

48.1 

7.6 

4.7 

-923.2 

-136.6 

Poitiers 

46.6 

0.3 

3.6 

-1517.2 

-186.7 

Schwarzenburg 

46.8 

7.3 

4.4 

-970.6 

-277.2 

Genova 

44.6 

9.0 

6.6 

-876.6 

-542.6 

Rome 

41.8 

12.5 

6.7 

-628.8 

-890.1 

Ebro 

40.8 

0.3 

5.0 

-1695.5 

-836.8 

Leningrad 

60.0 

30.7 

5.8 

603.8 

1172.9 

Uppsala 

59.8 

17.6 

3.8 

-136.2 

1102.7 

Gorki 

56.1 

44.3 

6.9 

1518.8 

953.7 

Moscow 

55.5 

37.3 

8.0 

1096.0 

751.4 

Miedzeszyn 

52.2 

21.2 

6.6 

81.8 

245.3 

Rostov 

47.2 

39.7 

9.6 

1499.0 

-120.3 

Graz 

47.1 

15.5 

6.0 

-341.1 

-312.8 

ArtiGcial  Point 

— 

— 

9.5 

1000.0 

-750.0 
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Table  2.  Iteration  process  for  producing  map  shown  in  figure  S. 


foF2 

MHz 

10  iter. 

Ai  d 

20  iter. 

Ai 

d 

40  iter 

Ai  d 

X 

km 

Y 

km 

1 

3.1 

-4.304 

0.034 

-4.345 

0.011 

-4.360 

0.001 

-1463.5 

1093.8 

2 

7.0 

0.226 

-0.108 

0.374 

-0.030 

0.423 

-0.004 

-426.5 

532.6 

3 

6.0 

-0.156 

-0.013 

-0.210 

0.027 

-0.258 

0.004 

-1013.5 

337.5 

4 

5.2 

-0.863 

-0.232 

-0.625 

-0.065 

-0.528 

-0.007 

-1435.4 

359.6 

5 

7.1 

0.141 

0.008 

0.115 

0.010 

0.098 

0.002 

-684.6 

224.4 

6 

8.1 

0.093 

0.117 

-0.014 

0.026 

-0.050 

0.003 

-386.0 

13.9 

7 

6.5 

0.663 

-0.083 

0.694 

0.000 

0.696 

-0.001 

-1106.3 

207.1 

8 

6.2 

-0.627 

0.091 

-0.772 

0.049 

-0.855 

0.007 

-1320.9 

-58.1 

9 

7.6 

0.486 

-0.045 

0.566 

-0.030 

0.618 

-0.005 

-923.2 

-136.6 

10 

7.0 

0.420 

-0.039 

0.457 

-0.014 

0.491 

-0.004 

-1517.2 

-186.7 

11 

7.2 

-0.666 

0.100 

-0.746 

0.018 

-0.776 

0.003 

-970.6 

-277.2 

12 

8.7 

0.707 

0.039 

0.699 

0.003 

0.706 

-0.001 

-876.6 

-542.6 

13 

9.2 

1.282 

0.023 

1.274 

0.002 

1.269 

0.001 

-628.8 

-890.1 

14 

7.9 

0.641 

-0.007 

0.641 

0.002 

0.636 

0.001 

-1695.5 

-836.8 

15 

7.3 

-0.376 

-0.001 

-0.377 

0.001 

-0.378 

0.000 

603.8 

1172.9 

16 

5.8 

-1.507 

0.000 

-1.524 

0.008 

-1.539 

0.001 

-136.2 

1102.7 

17 

8.8 

0.117 

0.105 

0.028 

0.019 

0.003 

0.001 

1518.8 

953.7 

18 

9.5 

1.223 

-0.007 

1.254 

-0.013 

1.277 

-0.002 

1096.0 

751.4 

19 

8.8 

1.415 

-0.030 

1.438 

-0.003 

1.437 

0.001 

81.8 

245.3 

20 

10.0 

2.423 

-0.041 

2.462 

-0.008 

2.470 

0.000 

1499.0 

-120.3 

21 

8.7 

0.425 

0.061 

0.411 

-0.007 

0.431 

-0.002 

-341.1 

-312.8 

av  =  7.414 
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Table  3.  Plasma  frequencies  anu  their  coordinates. 


FN 

X 

Y 

z 

-0.5 

-300.000 

0.0 

120.0 

1.1 

-300.000 

0.0 

180.0 

2.3 

-300.000 

0.0 

230.0 

4.5 

-300.000 

0.0 

380.0 

-0.5 

-180.000 

120.0 

115.0 

2.3 

-180.000 

120.0 

220.0 

5.5 

-180.000 

120.0 

360.0 

-0.5 

-60.000 

0.0 

105.0 

1.1 

-60.000 

0.0 

160.0 

2.5 

-60.000 

0.0 

210.0 

6.3 

-60.000 

0.0 

340.0 

-0.5 

60.000 

120.0 

93.0 

1.1 

60.000 

120.0 

120.0 

0.8 

60.000 

120.0 

135.0 

2.0 

60.000 

120.0 

160.0 

3.3 

60.000 

120.0 

210.0 

7.0 

60.000 

120.0 

330.0 

-0.5 

180.000 

0.0 

85.0 

2.2 

180.000 

0.0 

110.0 

1.5 

180.000 

0.0 

130.0 

4.1 

180.000 

0.0 

210.0 

7.5 

180.000 

0.0 

320.0 

-0.5 

180.000 

120.0 

89.0 

3.2 

180.000 

120.0 

115.0 

2.0 

180.000 

120.0 

132.0 

3.3 

180.000 

120.0 

155.0 

4.8 

180.000 

120.0 

205.0 

8.2 

180.000 

120.0 

310.0 

-0.5 

300.000 

120.0 

85.0 

3.5 

300.000 

120.0 

110.0 

2.8 

300.000 

120.0 

130.0 

3.8 

300.000 

120.0 

150.0 

5.5 

300.000 

120.0 

200.0 

8.9 

300.000 

120.0 

300.0 

Form  Approved 
0MB  No.  070^188 


REPORT  DOCUMENTATION  PAGE 

PubNcraponlngbuRlantortMtcoNMtlonof  Moniudlanhestlinatedloawerag*  1  tiourpar  ratpons*.  Inctuding  Via  ttma  tor  lavtowtoig  kotiucllont,  aaarctikig  aoiMlng  dM  fouroat,  gMhtiIng  and 
riialiittlntogViad»tanaadad.aitocoiTiptollngaitoiavla»dngViacoMae«oiiofkitofmaaori.SattoeetninaritoraoafdlngV>hbuidana8l>na»ofanyoViafa»paelorvitocoiacltoiiorwotina<toti.^^ 
suggaeVlonstorraduclngtMsbuidaatoWasMngtonHaadquaitaisSa(Vtoaa.Dliaclo(alatorliito(inatlonOpatttlonaandRapoi1s.  1215Janinan  Davis  Highway,  Sulla  1204.  Aittigton.VA  22202-4302. 
andlothaOlllcaol  Managatnam  and  Budget.  PaparworkRadudlonPiolact  (0704-0168),  Washington,  DC  20S03. _ 


1.  AGENCY  USE  ONLY  (LstWfitonID  aREPORTOATE 

March  1991 

3.  REPORT  TYPE  AND  DATES  COVERED 

Final 

4.TTn£ANOSUBTTn£ 

5.  FUNDING  NUMBERS 

FTITING  OF  DISCRETE  IRREGULARLY  SPACED  DATA  WITH 
DIFFERENTIABLE  FUNCTIONS 

Application  to  Ray  Tracing  in  the  Ionosphere 

PE:  0602435N 

WU:  DN888715 

e.  AUTHOR(S) 

A.  K.  Paul 

7.  PERFORMING  ORGANIZATUN  NAME(S)  AND  ADORESS(ES) 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

Naval  Ocean  Systems  Center 

San  Diego,  CA  92152-5000 

NOSCTR  1405 

0.  SPONSORING/MONITORINQ  AGENCY  NAME(S)  AND  AOORESS(ES) 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

In  house 

1 1 .  SUPPLEMENTARY  NOTES 

12a.  OISTRIBUnON/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. 

13.  ABSTRACT  (MantmuTi  200  words) 


This  report  presents  a  new  approach  to  general  ionospheric  mapping  problems.  The  Htting  of  discrete  irregularly 
spaced  data  with  differentiable  functions  is  a  superposition  of  local  function,  'ne  contribution  of  an  individual  local  function 
to  the  total  variation  decreases  as  the  distance  from  the  center  of  each  such  function  increases.  The  method  can  be  applied  to 
many  kinds  of  data  representation  problems,  ranging  from  interpolation  of  missing  or  irregularly  spaced  data  in  time  series 
analysis  to  mapping  of  empirical  regular  or  irregularly  spaced  multidimensional  data.  The  application  of  this  method  will  be 
demonstrated  for  two  ionospheric  problems,  mapping  of  standard  parameters  and  ray  tracing  in  an  empirical  ionosphere. 


14  SUBJECT  TERMS 

ionospheric  mapping 
ray  tracing 

17  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICAnON 

OF  REPORT  OFTHBPAQE  OFABSTOACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 


IS.  NUMBER  OF  PAGES 

27 

laPRCECooi 

20.  UMITATX3N  OF  ABSTRACT 

SAME  AS  REPORT 


NSN  754001-280-5500 


SMdMdkimase 


N8N7S4»01-2aM800  Sttndirt  fcrni  288 

UNCLASSIFIED 


INITIAL  DISTRIBUTION 


Code  0012 

Patent  Counsel 

(1) 

Code  0144 

R.  November 

(1) 

Code  54 

J.  H.  Richter 

(1) 

Code  542 

A.  K.  Paul 

(20) 

Code  542 

B.  Rose 

(1) 

Code  542 

D .  Sailors 

(1) 

Code  542 

R.  Sprague 

(1) 

Code  542 

J.  Ferguson 

(1) 

Code  952B 

J .  Puleo 

(1) 

Code  961 

Archive/ Stock 

(6) 

Code  964B 

Library 

(3) 

Defense  Technical  Information  Center 
Alexandria,  VA  22304-6145 

(4) 

Naval  Research  Laboratory 
Washington,  DC  20375-5000 

(3) 

Office  of  Naval  Research 

Arlington,  VA  22217-5000 

(1) 

Department  of  Commerce 
Washington,  DC  20230 

(1) 

NASA 

Greenbelt,  MD  20771 

(2) 

NOAA 

Boulder,  CO  80303 

(5) 

Hanscom  AFB,  MA  01731 

(3) 

U.S.  Army  CECOM 

Fort  Monmouth,  NJ  07703 

(1) 

Los  Alamos  National  Laboratories 

Los  Alamos,  NM  87545 

(1) 

Center  for  Atmospheric  Research 
Lowell,  MA  01854 

(1) 

Boston  University 

Boston,  MA  02215 

(1) 

University  of  Illinois 

Urbana,  IL  61801 

(2) 

Utah  State  University 

Logan,  UT  84322-4405 

(1) 

GRC ,  Inc . 

Albuquerque,  NM  87106 

(1) 

Science  Applications  International 
McLean,  VA  22102 

(1) 

Southwest  Research  Institute 
San  Antonio,  TX  78284 

(1) 

SRI  International 

Menlo  Park,  CA  94025 

(2) 

